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Abstract. The Lee-Friedrichs model has been very useful in the study of decay- 
scattering systems in the framework of complex quantum mechanics. Since it is exactly 
soluble, the analytic structure of the amplitudes can be explicitly studied. It is shown 
in this paper that a similar model, which is also exactly soluble, can be constructed in 
quaternionic quantum mechanics. The problem of the decay of an unstable system is 
treated here. The use of the Laplace transform, involving quaternion-valued analytic func- 
tions of a variable with values in a complex subalgebra of the quaternion algebra, makes 
the analytic properties of the solution apparent; some analysis is given of the dominating 
structure in the analytic continuation to the lower half plane. A study of the corresponding 
scattering system will be given in a succeeding paper. 
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1. Introduction. 

Adler 1 ' 2 has developed methods in time dependent perturbation theory for the treat- 
ment of scattering and decay problems in quaternionic quantum theory. Although the 
spectrum of the anti-self-adjoint generator of evolution is quaternion definedf as positive, 
the resulting formulas for the scattering matrix display both left and right hand cuts. It 
has recently been shown 3 that such an anti-self-adjoint operator on a quaternionic Hilbert 
space has, quite generally, a symmetric effective spectrum, and in case the quaternion 
defined spectrum is absolutely continuous in [0, oo), the effective spectrum is absolutely 
continuous in (—00,00). In this case, a symmetric conjugate operator exists. In case the 
anti-self adjoint operator is the generator of evolution, the conjugate operator can be in- 
terpreted as a "time operator" ; it was shown 3 that there is then no evident contradiction 
to the definition of a Lyapunov operator 4 . 

The Lee-Friedrichs model 5 has been very fruitful in the study of decay-scattering 
systems in the complex Hilbert space 6 . The generalized states (elements of a Gel'fand 
triple, or rigged Hilbert space) associated with the resonances of the scattering matrix 
and the reduced resolvent, for the case in the the unperturbed Hamiltonian is bounded 
from below 7 as well as for the case for which, as in a model for the Stark effect, the 
"unperturbed" Hamiltonian has spectrum —00 to 00 8 , have also been studied using this 
spectral model. For the study of the consequences of the symmetrical effective spectrum 3 of 
a quaternionic quantum theory, as well as for the construction of physical models, it would 
be useful to construct an exactly soluble model of Lee-Friedrichs type in this framework 
as well. It is the purpose of this paper to construct such a model. 

We first utilize (Section 2) the formulation of the decay problem suitable for applica- 
tion of time-dependent perturbation theory, following the methods given by Adler 1 ' 2 and 
show that a soluble model can be constructed. In Section 3 we formulate the problem 
in terms of the Laplace transform of the Schrodinger evolution. The resolvent formalism 
cannot be used directly, since there is no natural imaginary direction associated with an 
anti-self-adjoint operator on a quaternion Hilbert space. Choosing a direction arbitrarily, 
and hence a subalgebra C(l,i) of the quaternion algebra H, the Laplace transform (and 
its inverse) of the Schrodinger equation is well defined. In the complex Hilbert space, for 
a time-independent Hamiltonian, this procedure is equivalent to the resolvent formalism; 
for a time dependent generalization of the Lee-Friedrichs model, it has been shown that 
it provides a useful basis for a well-controlled perturbation theory 9 . For the quaternionic 
analog of the Lee-Friedrichs model, we find that the Laplace transformed amplitude has 
a structure analogous to that of the resolvent form, but contains a linear combination of 
a function analytic in the transform variable z G C(l, i) and a function of the negative of 
its complex conjugate — z* , corresponding to the contributions of the right- and left-hand 
cuts. The analytic continuation of the amplitude into the lower half plane is carried out 
in Section 4, and some results are given for the structure of the functions there that can 



f For an anti-self-adjoint operator S, the (generalized) eigenvalue equation is Sf = 
fe\,e G ImH, e 2 = —1, where H is the quaternion algebra spanned over the reals by 
1, k, ij = k and i 2 = j 2 = k 2 = —1. One can choose the quaternion phase 1 of / such 
that e = i and A > 0. We call this the quaternion defined spectrum. 
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dominate the inverse transform. In a succeeding paper, some aspects of formal scattering 
theory are developed, and it is shown that this model provides a soluble scattering problem 
as well; it therefore describes a decay-scattering system. A summary and conclusions are 
given in Section 5. 

2. Spectral Expansion 

In this section, we follow the methods of Adler 1 ' 2 for the construction of an effective 
time dependent perturbation theory; assuming, as for the Lee-Friedrichs model, that the 
potential has no continuum-continuum matrix elements, we obtain a soluble quaternionic 
model for a decay-scattering system. 

In the framework of the theory of Wigner and Weisskopf 10 , the component of the 
wave function evolving under the action of the full Hamiltonian which remains in the 
initial state corresponds to the probability amplitude A(t) for the initial state to persist 
(the survival amplitude 11 ). Let \ipo) be the initial state and H be the quaternion linearf 
operator which generates the time evolution of the system. Then, 

A(t) = (Vo|e-^|Vo>. (2-1) 

Let us now suppose that H has the decomposition 

H = H + V, (2.2) 

where V is a small operator which perturbs the eigenstates of Hq. If we consider the 
eigenstates of Hq as the quasistable states of the system, (2.2) generates the evolution of 
an unstable system. In particular, if \ipo) in (2.1) is an eigenstate of Hq, the evolution of 
A(t) reflects its instability. 

Let us assume that there is just one eigenstate of H , with eigenvalue E , so that the 
spectral expansion of \ip(t)) = e~ Ht \ip ) is 

poo 

m)) = Me- iEot C (t) + / dE\E)e- iEt C E {t), (2.3) 

Jo 

where 

#o|V>o> = |V>o)^o (2.4) 

and 

(E\H \f)=iE{E\f). (2.5) 

We take initial values to be 

C (0) = 1, C E (0)=0. (2.6) 

| There are operators on the quaternionic Hilbert space with weaker linearity properties, 
such as complex linear or real linear 12 . We restrict ourselves here to the quaternion linear 
case, so that H(fq) = (Hf)q for all vectors / in the domain of H. 
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The transition from the eigenstate of Hq to the continuum is interpreted as the decay 
of the unstable system. The corresponding amplitude is 

(E\m) = e- iEt C E (t). (2.7) 

The total transition probability is 

poo poo 

/ \(E\m)\ 2 dE= \C E (t)\ 2 dE 

Jo Jo \*-°) 

= l-\C (t)\ 2 

The Schrodinger equation for the evolution, 

||V(*)> = -H\m), (2-9) 



with (2.3), (2.4) and (2.5), becomes 

poo 

\i>o)e- lEot C' (t) + / dE \E)e- iEt C' E {t) 
Jo 



-V\iP )e- tEot C (t) - / dE V\E)e- iEt C E {t). 

Jo 



(2.10) 



We take the scalar product in turn with and (E\ to obtain the coupled differential 
equations 2 

C '(t) = -e^MVl^e-^Coit) 

- / dE e iEot {^ \V\E)e- iEt C E (t), 
Jo 



and 



C E '(t) = -e lEt (E\V\i;o)e- lEot Co(t) 

dE' e iEt {E\V\E')e- iE,t C E ,(t). 



i>00 

L 



I AEt I Tp\-Tr\T?'\ „-iE' 't n ■ (2-12) 



'0 

Integrating these equations over a very short interval in t, covering t = 0, one sees 
that if we take Co(t), C E (t) to be zero for t < 0, the conditions (2.6) can be realized in the 
form C o (0_|_) = 1,Ce(0 + ) = 0, by adding a term 5(t) to the right hand side 2 of (2.11). 

We now assume a condition on V parallel to that of the Lee-Friedrichs model in the 
complex quantum theory, i.e., that 

(E\V\E')=0. (2.13) 

The imaginary complex part of (ipo\V\ipo), as for the Lee-Friedrichs model of the complex 
quantum theory, contributes a level shift (i.e., it appears only in the combination Eq + 



v, where iv is the complex part of V^[c/.(2.16)]), but the imaginary quaternionic part 
contributes to the quaternionic phase of the (sub-asymptotic) decay amplitude*. We shall 
retain the term (ipQ\V\ipo) here to display its effect. The condition (2.13) is essential to the 
solubility of the model; it corresponds to neglecting final state interactions, often (as in 
radiative processes) a good approximation to the physical problem in the complex quantum 
theory. Eqs. (2.11), (2.12) then become 

C '(t) = -e^^olVl^e-^Coit) 

r°° (2 14) 

-/ dE e tEot (iJo\V\E)e- lEt C E (t) + 5(t) 
Jo 

and 

C E '(t) = -e lEt (E\V\ifj )e- lEot C (t). (2.15) 

To make use of the methods of complex analysis, we decompose all of the quaternion- valued 
functions in (2.14), (2.15) into their symplectic components 2 ' 12 ' 13 , i.e., 

Co(t) = C 0a (t)+jC (3(t) 
C E (t) = C Ea (t)+jC Ef3 (t) 

(E\V\ip ) = -tyo\V\E)* = V Ea +jV Ef ,, (2 ' 16) 
(ipo\V\ip ) = iv+ju, 

where, since (ipo\V\ipo)* = — (-^ol^lV'o)? v is rea l 5 an d u G C(l,z); the asterisk is used here 
to denote the quaternion conjugation (it then corresponds to complex conjugation on the 
complex subalgebra C(l,z)). Then, Eqs. (2.14) and (2.15) become 



CUt) = -ivC 0a (t)+u*e 2iEot C p(t) 

poo 

+ / dE e^ E °- E ^V Ea C Ea (t) 
Jo 

pOO 

+ / dEe^ Eo+E ^V E ^C E ^t) + 5(t), 
Jo 

C'opV) = -e- 2iEot uC 0a (t)+ivC p(t) 

poo 

- / dE e-^ +E ^V Ep C Ea (t) 
Jo 

poo 

+ / dE e-^ E °- E ^V Ea C Ef3 (t) 
Jo 



(2.17) 



and 



* Since we have already chosen the quaternionic phase of \ipo) so that the eigenvalue 
of Ho is iE, with E > 0, there is no further freedom to reray {ipo\V\ipo) so that it lies in 
C(l,i). 
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C' Ea {t) = -e^ E - E ^V Ea C 0a (t) + e^ E+E ^V^C p(t) 
C' E0 (t) = -e-« E+E ^V Ep C 0a (t) ~ e-^ E - Eo)t V^ a C om . 
We now define the Fourier transforms 2 



(2.18) 



C a(t) = ~— / dwe^-^Coaiu) 
2l™ J-oo 

1 f°° 

C Ea (t) = -— dwe< E -^C Eoi {u) (2.19) 
2m J_ 00 

C Ef3 (t) = — / due-^+^CE^u) 
2m J_ 00 

1 f°° 

iS(t) = — / dwe !tE »- w)t 

27TZ 

Note that it follows from the conditions Co(t), C E {t) = for t < 0, that the complex- 
valued functions Co a (u), Cop(u), C Ea (u>), C E p(uj) are analytic in the upper half complex 
w-plane (upper half plane analyticity is equivalent to the condition that Co(t), C E {t) vanish 
for t < 0). 

With the definitions (2.19), Eqs. (2.17) and (2.18) become 



and 



-i = i(oj Eq v)C 0a (uj) + u*C p(u) 

+ / dE(-V£ a C Ea (u>) + V^C Ep {u)) 
Jo 

= i(uj + E + v)C 0(3 (uj) + uC 

poo 

+ / dE(V Ef3 C Ea (cv) + V Ea C Ef3 (u)), 
Jo 

C Ea (uj) = -^—{V Ea C 0a (oj) + V Ef3 C p(uj)) 

Hi — UJ 

C ep (uj) = -!—{y E pC Qa {uj) - v* a C Q p{u)), 

Hi "T Ul 



(2.20) 



(2.21) 



where the division is well-defined for u in the upper half plane. Substituting (2.21) into 
(2.20), one obtains the closed form 

h{u)C 0a {u) + g{-uyC Qp {u) = 1 

g(uj)C Oa (u)-h(-uj*)*Cop(uj) = 0, [ - ' 

where 
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g(oj) = -iu+ / d£ V Ea V E( 3 ( -=^— + ) 
Jo \E - ou E + uj J 



(2.23) 



Note that h(ui) and g(ui) are analytic in the upper half plane; it therefore follows that 
h(— uj*)* and g(— uj*)* are also analytic functions of uj in the upper half plane. 

The Eqs. (2.22) provide a closed solution for the coefficients C 0a (uj), Cap(u) in the 
domain of analyticity Im uj > 0. The time dependent coefficients Co a (t),Cop(t) can then 
be obtained by inverting the Fourier transform, integrating on a line infmitesimally above 
the real axis in the uj plane (implicit in the relations (2.19)). We wish to emphasize that 
unlike the Lee-Friedrichs model with semi-bounded unperturbed Hamiltonian, the cut is 
not just a half line. In accordance with what we have found for the effective spectrum of 
an anti-self-adjoint operator in ref. 3, these results reflect the presence of both a left and 
right hand cut 14 . The integrals can be estimated by lowering the contour into the lower 
half plane (the second Riemann sheet) 8 , where the analytic continuation of the first sheet 
functions are well-defined, to obtain the contribution of leading singularities. We study 
this analytic continuation in Section 4. 

3. Laplace Transforms 

In this section we study the Laplace transform of the amplitude A(t) defined in (2.1). 
As we have pointed out above, the direct application of resolvent methods is complicated by 
the quaternion structure of the Hilbert space. However, the Laplace transform nevertheless 
brings the equations to a form which is similar to that of the resolvent method in the 
complex Hilbert space. 

The amplitude A(t) satisfies the equation 

-^A(t) = (fo\(H + V)\m) (3J) 
= iE A(t) + (V>o| W(*)>, 

where we have used (2.4). 

Inserting the completeness relation 

poo 

l = |Vo)(^o|+ / dE\E)(E\ 
Jo 

between V and \ip(t)}, this is 

- m A(t) = iEoA(t) + J Q dE{to\V\E)A E (t), 
+ (^o\VMA(t), 

where 
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A E (t) = (E\m)- (3-3) 

The amplitude A E (i) satisfies 

~A E (t) = {E\H + V\m) (3 _ 4) 

= iEA E (t) + (E\V\ii }A(t), 

where we have used the completeness relation again, taking into account the restriction 
(2.13) on the matrix elements that constitute the soluble model, i.e., that V connects the 
continuum only to the discrete eigenstate. We now define the Laplace transform: 



A(z) = [ dt e lzt A(t), (3.5) 
Jo 

and a similar expression for A E (z). Clearly, A(z) exists and is analytic in the upper half 

z-plane, since by the Schwarz inequality 2 \A(t)\ < \\ipo\\ \\e~ Ht ipo\\ = ll^oll 2 , i.e., A(t) is 
bounded. The upper half plane analyticity of A E (z) follows from that of A(z), as we shall 
see below. 

The inverse Laplace transform follows as for the usual complex function space, since 
the quaternion structure of A(t) does not affect the left transform of A{z). For C a line 
from — oo + ie to +oo + ie, above the real axis in the z-plane, we assert that 

A(t) = — [ dz e~ lzt A(z). (3.6) 
2tt Jc 

For t < (A(z) — > for Im^ — > zoo), the contour can be closed in the upper half- plane, 
so that, by upper half-plane analyticity, A(t) = 0. For t > 0, 



A(t) = — [ dze~ lzt r dt'e lzt ' A(t') 
2tt Jc Jo 

-i poo poo 

= _L / dt ' e «t-s) / dx e -i<t-t') A (t>), 

2tt Jq J-oo 



(3.7) 



and, since J*^^ dx e lx( - 1 1 = 2n5{t — t'), the assertion is verified. 
The Laplace transform of Eq. (3.2) isb] 

/>oo 

1 + izA(z) = iE A(z) - dE V£ a A E (z) 

Jo 

POO 

+ j dE V Ep A E (-z*) +ivA(z) +juA(-z*), 
Jo 

\ Note that the first term on the left hand side is due to integration by parts, i.e., 

/oo o poo 

dt e tzt (—^A(t)) = -e tzt A(t)\^ + iz J dt e tzt A(t) = 1 + izA(z) 



(3.8) 
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and of Eq. (3.4), 

izA E (z) = iEA E (z) + V Ea A(z) +jV E pA(-z*). (3.9) 

These formulas reduce to the results for the complex case if Ve/3 and u are zero. The 
contribution of (i^>o\V\ipo) is then, as is clear from (3.8), just an energy shift. 

In terms of the components A a (z), Ap(z), Aeoi{z) and A E a(z), where (we shall see 
below that the a components are analytic functions of z, and the (3 components, in spite 
of the notation, are analytic functions of z*) 



A{z) =A a (z) +jAp(z), 
A E (z) = A Ea (z) + jA E p(z), 



(3.10) 



one obtains 



pQO 

1 + izA a (z) = iE A a (z) - dE V Ea A Ea (z) 

Jo 

- dE V ES3 A ES3 (-z*) + ivA a {z) - u*Ap(-z*), 
Jo 

/•OO 

-iz*Ap{z) = -iE A[3(z) - / dE V Ea A E/3 (z) 

Jo 

poo 

+ dE V Ef3 A Ea (-z*) - ivAp(z) + uA a (-z*), 
Jo 

and 

izA Ea (z) = iEA Ea (z) + V Ea A a (z) - V Ep A p (-z*) 
-iz*A Ep (z) = —iEA E p(z) + V Ea A p {z) + V Ep A a {-z*) 



(3.11) 



(3.12) 



We now study the relation between the amplitudes A(z), A E (z) and Co(z), C E (z), 
and show that the equations (3.8) and (3.9) are equivalent to (2.20) and (2.21). 
The transforms defined in (2.19) have the property 

C Ea (t) + jC Ef3 (t) = -—e lEt / dw{e- wt CE a {w) + e Mt 3C E ft{w)}, (3.13) 

2m J_ 00 

a relation valid for the CE (t) coefficients as well. The transforms act with opposite sign 
of the frequency on the a and (3 parts of the C-amplitudes; they are, however, analytic 
in the upper half u;-plane by construction. On the other hand, the transform (3.5), also 
analytic in the upper half z-plane, for 

A(t)=A a (t)+jA l3 (t), (3.14) 

has the form 

p OO p OO 

A{z)= I dt e tzt A a (t)+j dte- izH Af3{t), (3.15) 
Jo Jo 
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and hence, with (3.10), we see that 



/•oo 

A a (z) = / dt e izt A a (t) 
Jo 

/•OO 

A p (z) = / dt e- iz * l A^t). 
Jo 



(3.16) 



The second of (3.16) is an analytic function of z* in the upper half z-plane (for which 
—z* is also in the upper half plane). 

This form is, in fact a general one 15 for quaternion-valued left analytic functions of 
a variable in the complex subalgebra C(l,i). To see this, consider the Cauchy-Riemann 
relations for the existence of a left-acting derivative with respect to a variable in the 
complex subalgebra for 



f(z) = f a (z)+jfp(z) 

= u a + iv a +j(up + ivp). 



df du a .dv a .( dup . dvp 



Then, 



dx dx dx \ dx dx 
1 df_ _ Idu^ dva_ _ . ( 1 Sup dvp \ 
i dy i dy dy \i dy dy J 

so that if the left derivative of f(z) exists, 

du a dv a dv a du a 



dx dy dx dy 

and 

dup _ _dvf3_ dVf3_ _ dup 

dx dy dx dy 



(3.17) 



(3.18) 



The relations (3.18) correspond to anti-Cauchy-Riemann conditions, implying that fp 
is an analytic function of z* [note that with the convention we have used, i.e, of writing j to 
the left of the (3 part of the function, right analyticity results in normal Cauchy-Riemann 
conditions for both parts; later we shall study a function with j appearing to the right of 
the /9-part, for which the normal Cauchy-Riemann conditions for left analyticity apply to 
both parts] . 

From (2.3), we have the relations 



A{t)=e-^ t C {t), 
A E (t) = e- iEt C E (t). 
It then follows from (2.19) that the Laplace transform of A(t) is 
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j^^fJ^-j C °f\). (3.20) 

The functions Co a (u), Cop(u>) are analytic in the upper half w-plane, and for smooth 
functions Co a (t), Cop(t), vanish for Im u — > oo. Hence, one can complete the contour in 
(3.20) to obtain 

A a (z) = iCn a (z) 

(3.21) 

A (z) = -iC Qf) (-z*), 

consistent with the upper half plane left analyticity of A(z). The functions A Ea (z), Aep(z) 
are related in the same way to Ce<x(z), C E p(— z*). 

Substituting the relations (3.21) into (3.11) and (3.12), one obtains 

-i = -i( z -E - v)C 0a (z) + u*C Gf3 {z) 

/■OO 

+ / dE(-V Ea CEa(z) + V Ef3 C E p(z)) 
Jo 

= i(E + v- z*)C 0f 3(-z*) + uCoa(-z*) 

/■OO 

+ / dE{V E0 C Ea {-z*) + V Ea C E(3 {-z*)), 
Jo 



(3.22) 



and 



C Ea (z) = {V Ea C 0a (z) + V^C 0f3 (z)) 

Hi — Z 

CE f3 (-z*) = ¥ ^(VEpCo a (-z*)-V^Co (i (-z*)). 



(3.23) 



The replacement of — z* by z in the second of each of (3.22) and (3.23) corresponds to an 
analytic continuation of these equations inside the upper half plane, where the functions 
are analytic. This replacement brings these expressions into correspondence with (2.20) 
and (2.21). 

Let us now return to the quaternionic equations (3.8) and (3.9). Substituting the 
second into the first to eliminate the amplitude A E (z), we obtain 

1 + ih(z)A(z) = -ijg(-z*)A(-z*), (3.24) 

where h and g are the functions defined in (2.23). Note that although it appears that 
g(z) , as defined by (2.23), is an even function of its argument, it is , in fact, defined in 
the upper half plane, and has a cut along the real axis, so that the formal reflection of 
z takes the function out of its original domain of definition. The domain of definition of 
the function can be extended to the lower half plane by analytic continuaton, which we 
shall discuss below. The same remark is applicable to the function h(z), which appears, in 
(2.23), to have a trivial reality property, i.e., that h(z)* = h(z*); this operation, however, 
also shifts the variable from the upper to the lower half plane. The value of the extension 
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of this function to the lower half plane also involves analytic continuation, which we shall 
discuss below. 

If Ve/3 vanishes, then g(z) vanishes, and the result (3.24) coincides precisely with 
the well-known result for complex quantum theory, usually obtained by means of the 
method of resolvents 6 ' 7 . The quaternionic part of the transition matrix element couples 
the amplitude to its reflection in the imaginary axis of the complex energy plane, and 
corresponds to the effect of the negative energy part of the effective spectrum of the anti- 
self-adjoint Hamiltonian 5 . 

Since h(z), g(z) and A(z) are analytic in the upper half plane, Eq.(3.24) can be shifted 
to the point —z*, providing the second coupled equation 

1 + ih(-z*)A(-z*) = -ijg(z)A(z). (3.25) 

Solving for A(— z*) from the (3.25), and substituting into (3.24), one obtains the relation 

H(z)A(z) = iG(z), (3.26) 

where 

H(z) = h(z)h(z)+g(z)g(z) (3.27) 

and 

G(z) = h(z)+g(z)j, (3.28) 

and we have used the notation 

h(z) = h{-2?)\ g{z)=g{-z*Y. (3.29) 

The functions h(z) and g(z) are analytic functions of z if h(z) and g(z) are; hence H(z) and 
G(z) are analytic functions of z as well. The function g(z) contains the complex quantity 
u, the quaternionic part of (tpo\V\ij;o), which therefore contributes to the quaternionic part 
of the decay amplitude A{t). 

4. Analytic Continuation 

The inverse of the Laplace transform (3.6) can be well-approximated when there are 
singularities in the analytic continuation of the amplitude A(z) into the lower half plane, 
close to the real axis, by deforming the countour of the integration (3.6) below the real 
axis. Since, in Eq. (3.26), G(z) depends on functions which can be taken to have regular 
analytic continuation to the lower half plane, such singularities arise, in general, from zeros 
of H(z). 

We shall study in this section the analytic properties of h(z) and g(z) to see how zeros 
of H{z) could arise. To do this, it is convenient to put these functions into a somewhat 
simpler form by defining numerator functions in the integrals that can be continued from 
the interval [0, oo) to (— oo, oo). As pointed out in ref. 5, the effective negative spectrum 
of Hq can be represented by the identification (we use a double bar on the Dirac kets to 
indicate the spectral representation for which the argument may be positive or negative), 
for E > 0, 

(-E\\ = -j{E\. (4-1) 
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The matrix element 

V E = (E\V\il> ) (4.2) 

then has the extension 

V-e = (-E\\V\iPo) = -j(E\V\^o), (4.3) 

so that, for E > 0, 

V-Ea = Ve/3 

and 

V-ep = ~V Ea . (4.4) 

Since the limit of V-Ea as E — > + is not necessarily equal to the limit of Veo. as E — > 0+ 
(we do not require continuity across zero), (4.4) does not place any restriction on the 
relation between Veo. and Ve/3- 

With these identifications, it follows that 



\{E\\V\^ )a\ 2 f° ,„ \(-E\\V\ik) 12 



a 



oo E - z E + z 

oo I T /_ 1 2 







(4.5) 



so that we can write 

h(z)=z-E -v+ [ dE (4.6) 

J-oo E-z 

where 

X(E) = |V^ Q | 2 , (4.7) 

for E E (— oo, oo). 

In a similar way, since (by (4.4)) 

V-EaV-E/3 = —VEpVEa, 

we may write 

g{z) = dE - iu, (4.8) 

where 

F(£) = -Y(-E) = V Ea V E p. (4.9) 
Taking the imaginary part of (4.6), one obtains 



( 1 + f 

\ J — OO 



\mh(z) = fin.: I I + / r/£ ^ E \ 2 )• ( 1-10) 



so that there are no zeros of h(z) for z in the upper half plane (this is therefore true of 
h(z) as well). 



13 



Taking z onto the real axis from above, one finds, from (2.23) with z in place of u, 



h(X) 



X - E Q - v + n i 1 V Xa | 2 + P I Q °° dE ( %f - ) , for z - A > ; 

-X-E -v + m\Vxp\ 2 + Pj^dE ^^-\^efj, for z -> -A < 0, 



(4-11) 

where the principal part restriction is not necessary in one term of each of the integrals. 
Since |F_Aa| 2 = |^A/3| 2 5 this result follows, equivalently, from (4.6) as well. If we now 
assume that we have chosen a model (specified by v, u and the functions Veo. an d Ve&) for 
which X(E) is the boundary value of a function X(z) analytic in the lower half plane in 
a large enough domain below the real axis to cover the analytic structure we must study, 
h(z) has the analytic continuation in the lower half plane 

/oo X ( E) 
dE — ^ + 2niX(z). (4.12) 
-oo E-z 

We see that taking the limit of h 11 (z) onto the real axis from below, one obtains exactly 
(4.11); a term — 7ri|V\ a | 2 for z — > A > 0, or — 7rz|Vx^| 2 for z — > —A < 0, emerges from 
the integral, cancelling half of the additional term 2niX(z). Hence h(z) and h II (z) are 
functions analytic, respectively, in some domain in the upper half plane and in some domain 
of the lower half plane, with a common segment of the real axis on which they approach 
the same values; a simple application of the Cauchy theorem demonstrates that they form 
one analytic function, and hence are unique analytic continuations of each other. It then 
follows that 



bn/r"(>) = 2nRcX(z) + hn: ( ! + / dE Jr^L )• < I - 13) 



which, if KeX(z) > 0, to be expected if z is close to the real axis, can vanish for Imz 
negative. We shall assume for the remainder of the analysis, for simplicity, that X(E) <^ E 
for every E, so that the assumption that the root of h n {z) has small imaginary part is 
justified. 

In a similar way, one finds the analytic continuation of g(z) to be 

g n (z) = [ dE + 2niY(z) - iu, (4.14) 

J-oo E-z 

where we have assumed Y(E) to be the boundary value of an analytic function Y(z), 
analytic in the same domain as X(z). There are no model independent statements we can 
make about g II {z) since Y(E) has no definite sign or reality properties; we shall therefore 
further assume, for the purposes of our present discussion, that Veq Veo., an d that u 
is also small, so that the principal contribution to the contour integral (3.6) comes from 
values of z close to roots of the analytic continuation of the first term, h(z)h(z), in the 
expression (3.27) for H(z). With this assumption, only a small adjustment in the value z 
of a root of this term is required to cancel the contribution of the second term as well. We 
therefore study the zeros of h(z)h(z). 
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It follows from (4.12) that the root ( must satisfy 

X(E) 



Reh n (C) = = ReC( 1 - / dE 



\ J —oo 



oo 
oc 



l^-CI 



( 4 - 15 ) 

-E -v + j dE E lE _^ 2 ~ ^lmX(C), 

so that a zero of the real part can only occur at 

u = ReC = E + v. (4.16) 

If \v\ < Eo, it follows from (4.12) and (4.13) (and the assumption that X(E) is small) that 
a zero can only form in h 11 (z) below the positive real axis. The vanishing of the imaginary 
part then requires, according to (4.13), that 

ImC = -nReX(C) - -n\V^ a \ 2 , (4.17) 

where we understand that ujq = Eq + v is positive, and we have approximated X(Q by 
its value on the real axis. This result is obtained from (4.13) by noting that the integral, 
for small values of Imz, contains an approximation to a representation of the 5-function, 
i.e. (for e > 0), 

lim 7T = tt6(x). 

Hence, 

I Tm y I 

^ ty5(E- Rez); 



\E-z\ 

substituting this back into (4.13), with Imh 11 (z) = 0, one obtains (4.17). Using the same 
approximation, this result may be obtained directly from (4.12); the integral, for small 
negative Imz contributes an imaginary term which is approximately —niX(uJo). 

We now show that there is a second leading contribution to the inverse transform (3.6) 
when the contour is deformed into the lower half plane. The function 

h(z) = h(-z*)* = -z - E - v + [ dE (4.18) 



oc 



E + z 



is analytic in the upper half plane and (following the argument leading to (4.10) for h(z)) 
has no zeros for Imz ^ 0; its analytic continuation to the lower half z-plane, obtained 
following the procedure resulting in (4.12), is 

h U ( z ) = - z -E -v+ [ dE -2mX(z), (4.19) 

J-oo E + z 

where X(z) is a function analytic in the lower half plane, in a sufficiently large domain, 
with boundary value X(—X) on the real axis, for z — > A. Since, from (4.19), we have the 
relation 

Imh H (z) = -lmz(l+ [ dE ,f"^, 2 > ) -27rReX(z), (4.20) 
V i-oo \E + z\ J 
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h 11 (z) can vanish for Imz < 0. 
The real part of h (z) is 



Reh U (z) = -Rez( 1 - / dE *^ E \A -Eq-v 

dE E . L +2nlmX(z), 
-oo \E + z\ z 

and can therefore vanish at the root £ for 

ReC = -u = -(E + v). (4.22) 

From (4.20), using the procedure described above to obtain (4.17), we find that 

ImC= -tvX(ujo), (4.23) 

where we have approximated X(Q) by its value on the real axis, X(—u>o) = X(uo). As 
for (4.17), this result follows from (4.20) by using the approximation 

-Imz I dE *^ E \ 2 = ivX(-Rez) (4.24) 
J-oo \E + z\ 

for Im^ small. Hence the decay rate associated with the residue of the second pole is the 
same as that of the first. 

We wish now to estimate the contour integral (3.6) when the contour is deformed into 
the lower half plane ft- The integral is given by 

— / dze- lzt A(z)^ e ~ lZt —i GlI ( z ) + background, (4.25) 

where we neglect, for simplicity, the second term of H 11 (z) in the denominator, and the 
background contribution consists of the contour integral in the lower half plane below 
the pole contributions; with a suitable choice of the weight functions, for example, for 
X(z),Y(z) Hardy class in the lower half plane on lines with sufficiently large negative 
imaginary part, the background contribution can be made as small as we wish (as in the 
Pietenpol model 8 in the complex Hilbert space), given e > 0, for t > e. 



(j For large real part, the denominator function H(z) goes as 0(z 2 ), and the numerator 
G(z) as O(z). At the ends of the long rectangle that brings the integral to the lower half 
plane, the contribution of the vertical contour is therefore bounded by x~ x J dy e yt , where 
y is the (negative) imaginary part of z, x is the very large real part of z, and the integral is 
taken between zero and the distance into the lower half plane that the contour is carried. 
These contributions therefore vanish for \x\ — > oo, even for t — > 0. 
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To do this, we must also study the analytic continuation of 

g(z) ee g(-z*y = dE + iu*. (4.26) 

J-oo E + z 

Using the property Y(—E) = —Y(E) : it is convenient to rewrite (4.26) as 

g(z) = f°° dE^± + iu*. (4.27) 

J-oo E-z 

The analytic continuation is then given, following the procedure used to obtain (4.12) and 
(4.14), by 

dE — + iu* + 2niY(z), (4.28) 
-oo E-z 

where Y(z) is a function analytic in the lower half-plane with the same domain of analyt- 
icity as X{z), Y(z), and with boundary value Y*(X) in the limit z — > A onto the real axis. 
At the zero ( of h IT (z), close to the positive real axis, near E = uj , 

/°° Y*(E) 
dE - K ' + iu* + iviY* (u ) . (4.29) 
-oo E - UJ 

At the root close to the negative real axis, near E = — ujq, 

g n (C) *P [°° dE + iu* + 7vtY*(-iv ) 



' — OO 

oo 



OC 



dE ^ ( E > + iu* - TviY*(u ). 
E - ujq 



(4.30) 



We are now in a position to evaluate the pole contributions to (4.25). The residue at 
C contributes a term which is approximately 

^V'(C) = e-«* f 1 + , (4.31) 

where we have neglected small terms in the residue of h II {z). 

The zero of the denominator in the neighborhood of z ~ Q is due to the factor h (z) 
in H II {z). In this case, however, the first term of G II (z) vanishes, so that the contribution 
of the negative energy pole is 

_ e iCt9_JQ.j (4.32) 
h n (C) 

the quaternionic reflection of the resonance below the positive part of the real line. The 
exponential decay rate of both of these contributions is given by (4.17), i.e., for 

C = «b-*(7/2), (4.33) 
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the decay rate 7 is given by 

7 = 27r|K, oa | 2 . (4.34) 



The oscillation frequencies of the two terms have, however, opposite sign. 

From (4.12) and (4.19), we see that h H (C) and h (C) are both well-approximated by 
-2u ; we write (4.29) and (4.30) as 

/'(C)=& + 7uy> ) (4.35) 

and 

g H {Q = b-mY*{u ), (4.36) 

where 

/°° Y*(E) 

dE y ' +iu*. (4.37) 
-00 E-uj 

In terms of these quantities, 

Ait) = e" lCt - ^—\e-^\b + 7riY*(iv )) - e**\b - 7riY*(u ))Xj, (4.38) 

2uj I J 



and hence 



\A(t)\ 2 ^e-^\l + ^[\b\ 2 + n 2 \Y(u )\ 2 

1 (4.39) 

- (\b\ 2 -tv 2 \Y(uj )\ 2 ) cos2uj t - 2Re(Y(uj )nb sm2uj t] i; 



the term in square brackets is non-negative. It is clear from (4.38) that in the limit t — > 0, 
the approximate expression for A(t) does not approach unity (the approximations we have 
made in arriving at this expression do not affect the quaternionic amplitude in this order) . 
We see, moreover, that |A(t)| 2 approaches a value slightly larger than unity as t — > 0, in 
agreement with the usual result 6 , implied by the fact that the exact initial decay curve 
starts as 1 — 0(t 2 ) in the neighborhood of t = (this follows from the anti-self-adjoint 
property of H). We have argued, however, that since the background contribution can 
be made as small as we wish, the pole contributions must be as accurate as we wish, for 
t > e > 0. The non- cancellation of the quaternionic pole residues is an indication that the 
series expansion of A(t) for very small t is not smoothly connected to the non-perturbative 
exponential functions which can represent A(t) with, in principle, arbitrary accuracy for 
non-vanishing t. We shall discuss this phenomenon a little further in the next section. 

The probability of survival of the unstable state has small oscillatory interference terms 
of twice the frequency (Eo + v), with coefficients proportional to the square of VeoVe(3'i 
the oscillation in the survival probability could conceivably provide an observable residual 
quaternionic effect. Note that these oscillations do not correspond to those of the time 
dependent perturbation expansion for small t which lead to convergence to the Golden 
Rule 2 . 
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If the complex part of (V'ol^lV'o) is large and negative , i.e., v < —E , the corre- 
sponding complex Hilbert space problem develops a bound state due to the interaction 
(the complex pole moves to the left and onto the negative real axis). For the quaternionic 
problem, however, the pole condition now becomes 

ui = -ReC- = E + v. (4.40) 

The pole is, in this case, under the negative real axis, at z = where (with the same 
argument leading to (4.17)) 

ImC- = -*\V Ul p\ 2 . (4.41) 

This very small width corresponds to a resonance that is an almost-bound state. The 
occurrence of a resonance in quaternionic quantum theory replacing the bound state, in a 
corresponding situation, in the complex quantum theory, has been discussed by Adler 1 ' 2 
both for the one-dimensional and spherically symmetric three dimensional potential prob- 
lems. 

The residue of this pole at — u\ contributes to the survival amplitude a term given 
approximately by 

e-«-«(l + I^A (4-42) 

where h (£_) = 2u>\. As in the previous case, there is a second pole, now under the 
positive real axis, for which — z* = i.e., z = — C* = C-- The contribution of this pole 
to the survival amplitude is approximately given by 

_ e iCltt%Aj, (4.43) 

h"(C-Y 

corresponding to a resonance with positive frequency and very small amplitude, which is 
the quaternionic reflection of the almost-bound state on the negative part of the real axis. 



5. Conclusions 

We have shown that a soluble model for the description of a quaternionic quantum 
mechanical system undergoing decay according the description of Wigner and Weisskopf 10 
can be constructed in the same way as the Lee-Friedrichs model 5 of the complex quantum 
theory. 

The existence of the effective spectrum 3 in (— oo, oo) of the anti-self-adjoint Hamilto- 
nian implies that the problem of the analysis of the complex structure of the amplitude is 
quite similar to that of the Pietenpol model 8 , an idealization of the Stark Hamiltonian in 
complex quantum theory. In this model, the spectrum of the "unperturbed" Hamiltonian 
has an absolutely continuous part in (— oo, oo), and the Wigner- Weisskopf theory achieves 
an almost irreversible law of evolution (there is no branch point to contribute to long time 
deviations from exponential behavior). 

Carrying out an approximate analysis (assuming the discrete-continuous matrix ele- 
ments of the the perturbing potential are small, and that those of the quaternionic part 
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are even smaller) of the analytic properties of the (analytically continued) amplitudes in 
the lower half plane, we find that for perturbations resulting in an energy shift of the real 
part that does not become negative, the unperturbed bound state becomes a resonance 
of the usual type (although the amplitude has a quaternionic phase); there is, however, a 
reflection resonance of small amplitude, with negative frequency, below the negative part 
of the real axis. This reflected resonance results in an oscillatory interference term which 
could provide an observable quaternionic residual effect. If, on the other hand, the energy 
shift induced by the diagonal part of the perturbation moves the pole to a position below 
the negative part of the real axis, in place of the bound state that would be found in the 
corresponding complex theory, one finds a rather narrow resonance which becomes a bound 
state in the limit in which the quaternionic part of the potential vanishes. Its quaternionic 
reflection on the positive part of the real axis is a positive frequency resonance of small 
amplitude which vanishes, along with the width of the negative energy resonance, with 
the quaternionic part of the potential. The occurrence, in quaternionic quantum theory, 
of a resonance replacing a corresponding bound state in the complex quantum theory, has 
been discussed by Adler 1 ' 2 both for the one-dimensional and spherically symmetric three 
dimensional potential problems. 

As we have seen in the previous section, the quaternionic part of the pole residues at 
C and ( = —(* do not cancel to leave an approximate unity residue as t — > as one might 
expect. The non-cancelling part comes from the behavior of the quaternionic part g(z) 
of the amplitude near the real axis, in particular, from a numerator function (absorptive 
part) which emerges from the integral due to the 5-function that forms as z approaches 
the real axis. 

To understand this behavior, consider the Pietenpol model of the complex theory 8 . 
For a Hamiltonian of the form H = H + V, where H has continuous spectrum from 
— oo to oo, and a bound state cf> embedded at E > 0, the singularities of the function 
{<f>\R(z)\<f>) , where R(z) is the resolvent (z — H)~ x , are determined by 

M*)<«*)|0) = 1, (5.1) 

where 

,oo ,/ 

h(z) = 7, - Rn - V + I — 

E-z 



!,( = ) = z-E -v + / ' Hf, (5.2) 



— oo 



and v is the expectation value of V in the state 0. The analytic continuation of this 
equation to the lower half plane, as for eq.(4.12), has the form 



h n (z)=z-E -v+ f dE + 2niX(z), (5.3) 

J-oc E ~ z 

where X(E) = \{E\V\(j))\ 2 , the boundary value of the analytic function X(z) . Following 
the arguments we have given above, one sees that h 11 (z) = for 

z = C = -7riX(u ), (5.4) 
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for u)o = Eq + v. This result corresponds to a nonperturbative width for the decaying 
system. On the other hand, the time dependent perturbation expansion of the amplitude 
A(t) results in 

e^ Eo+v)t A(t) = 1 - I(t) + . . ., (5.5) 



where 



1(f) = dE X i E) U{E -E -v)t + l- e *(^-)*}. 
It then follows that 

= 1 - 27rX(£ + u)* 



(5.6) 



(5.7) 



providing a perturbative width in agreement with the nonperturbative width obtained in 
(5.4); the condition for this result, corresponding to the Golden Rule, is that t be greater 
than 1/Eq. However, as we have pointed out above, the non-perturbative calculation of 
the time dependence of the amplitude, as represented by the pole residues, can be made 
as accurate as we wish. The perturbation expansion, in low order, is therefore not an 
accurate representation of the decay law, in this case, for any t > e > 0. 

A similar phemonenon occurs for the spectrum of the Stark effect. The bound states 
of an atomic system, such as hydrogen, are destroyed by the Stark potential (see ref. 8 
for a discussion) for any electric field, no matter how small. The perturbative treatment 
of the energy levels, however, results in a very good approximation of the level shift from 
the initially bound states to the peaks of the spectral enhancements in the continuous 
spectrum 16 . 

Some development of a formal scattering theory and the properties of this model for 
the description of scattering will be given in a succeeding publication. 
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